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Abstract 

CO ■ The pentaquark width is calculated in QCD sum rules. The higher dimension 

operators contribution is accounted. It is shown, that Tq should be very small, 
' less than lMev. 

' 

O ■ PACS: 12.39 Dc, 12.39-x, 12.38 

^: 

Ph. The status of + , predicted in 1997 by D.Diakonov, V.Petrov and M.Polyakov pQ 

in the Chiral Soliton Model, till now is doubtful. Few years ago narrow exotic baryon 
resonance + with quark content @ + = uudds and mass 1.54 GeV had been discovered 
by two groups [3 El- But the question is open until now, during last two years some 
groups confirm the pentaquark + existence, while other see null signal. Moreover, last 
year some groups, which have seen pentaquark, in the new experiments with higher 
statistics reported null result for pentaquark signal (CLAS experiments on hydrogen 
and deuterium jl], BELLE j^j) but at the same time DIANA and also LEPS, SVD- 
2 confirm their results with higher statistic(see [7j for the review). (Some theoretical 
reviews are given in IH] ) • 

So yet one can say only that if pentaquark exist, it should be a narrow state. 
Experimentally, only an upper limit was found, the stringer bound was presented in 
[3]: T < 9MeV. The phase analysis of KN scattering results in the even stronger 
limit on V JUj, T < lMeV. A close to the latter limitation was found in [llj from the 
analysis of Kd —>■ ppK reaction and in ^2] from K + Xe collisions data [3 . Also [Sj 
from the negative result of the experiment give the upper limit for pentaquark width 
less than QAOKeV 

In the paper ^3], [H| it was shown, that if pentaquark is genuine state it width 
should be strongly parametrically suppressed. It is necessary to note that this is 
general statement and does not depend of the choice of the pentaquark current (without 
derivatives), but at the same time it is significantly based on the assumption, that the 
size of the pentaquark is not larger, than usual hadronic. The main goal of this paper is 
to find numerical estimation of pentaquark width by use the method offered in |13j.|14j. 



X 



1 



Part 1. In the papers ^3], [2j it was shown, that pentaquark width should be 
suppressed as r ~ Ci 2 (0\qq\0) 2 , (for any current without derivatives). Later, in the 
short paper the first non-vanishing operator (dimension d — 3) contribution was 
calculated and the sum rule was considered numerically. It was shown that pentaquark 
width is suppressed numerically also and the width of the pentaquark width was esti- 
mated to be less 1 Mev. In this paper we will discuss this sum rules in more detail and 
also the contribution of the operators of the higher dimensions will be accounted. Let 
us shortly remind the main points of the method. We start from 3-point correlator 

n„ = / e^*-^<0 I Ve(x)jl(y)U0) | 0) (1) 

where rj n (x) is the neutron quark current (TB], (r] n = e abc (d a Cj ll d b ) r y 5 'y^u c ), 

(0 | r] n | n) = X n v n , (v n is the nucleon spinor), t]q is an arbitrary pentaquark current 

(0 | r/0 | 6 + ) = XgVg and = s^^u is the strange axial current. 

As an example of rjg one can use the following one (see , where it was first offered, 

and also [TH], where the sum 2-point rule analysis for this current was discussed): 

J A = e abc e def e 9cf {u aT Cd b ){u dT CYl5d e )Ycs g (2) 

and we will use it farther to obtain numerical results. (Of course there is large 
number of the another currents, for example see [T^], where 2-point correlators was 
analyzed very carefully, taking into account operators up to dimension 13 and direct 
instanton contribution). 

As usual in QCD sum rule the physical representation of correlator (1) can be 
saturated by lower resonance states plus continuum (both in t]q and nucleon channel) 

K hVS = (° I Ve I I U I n )( n \Vn I 0) 2 1 2 2 1 n + cont. (3) 

where P2 = Pi — q is nucleon momentum, m and mg are nucleon and pentaquark 
masses. 

Obviously, in the limit of massless kaon 

(0 + I j, I n) = g£ n v n [g^ - Yl^ e (4) 

where axial transition constant gf n is just we are interesting in (the width is pro- 
portional to the square of this value). Such a method for calculation the width in QCD 
sum rules is not new, see, e.g. j2D]- in the case of massive kaon the only change is in 
denominator of second term in r.h.s of the eq. (4), i.e. 

(0 + I 3n I n> = gtv n IgT - -f^n) ^ ^ (5) 
V 1 ~ m k/ 

It is clear that the second term vanishes at small q 2 . 

Substituting (0 | r\ n \ n) = X n v n , and (0 | rje \ 6 + ) = XgVg in eq (3) and take the 
sum on polarization one can easily see, that (in the limit of small q 2 ) correlator (1) is 
proportional to gf n . 
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nf V s = XnXggtn- 5 j(-2ftrf 7 5 + ••••) (6) 

where dots in r.h.s mean other kinematic structures (proportional to q e.t.c). From 
(6) one can easily obtain sum rule for axial constant gf n . For our sum rules we will 
use invariant amplitude just at the kinematical structure because, as it was 

discussed in [21], [22], [23] the choice of the kinematic structures with maximal number 
of momentum lead to better sum rules. So we obtain the following sum rules 

A n A^ n e-«/ M »+™«) = (-1/2)5^11^ (7) 

where Bg, B n mean Borel transformation on pentaquark and nucleon momenta cor- 
respondingly, and continuum extraction is supposed. 

By use of the equation of motions the eq. (4) close to the mass shell can be rewritten 

(0 + I jfi I n) = gtv n (V + ^±^V) l5 v e (8) 

At the same time, the second term in (4,8) correspond to the kaon contribution to 
9 — n transition with lagrangian density L = igg n kVnl 5 'V8 < f ) k, so one can write 

{° + I fa I n ) = Qenk 2 ^ 2 V n ^Vg (9) 

Comparing (9) and (8) one can found (if we for a moment neglect the kaon mass) 

genkfk = (m n + mg)g£ n (10) 

This is the analog of the Golderberger-Trieman relation. Of course the accuracy of 
this relation is about the scale of SU(3) violation but as estimation of the value of gg n k 
it is enough good. In ^1] some general properties of correlator (10) and corre- 
spondingly sum rules for gg n was obtained. First of all, it was shown, that correlator 
(1) vanishes in the chiral limit for any pentaquark current without derivatives. That 
means, that first non-vanishing contribution to sum rules give the operator with di- 
mension d = 3 (quark condensate), so axial constant gf n should be proportional to the 
quark condensate. An examples of corresponding diagrams are shown on the Fig.la,b. 
But as was also shown in these papers, diagrams like those on fig. la (i.e. without hard 
gluon exchange) can not contribute to the sum rule. The reason is that such diagrams, 
as one can easily check, are expressed in terms of the following integrals 

e i(pix- m ) axa y = / e e d ± xdH (n) 



((x - y) 2 ) n (x 2 ) m J (x 2 ) m (t 2 ) n 

It is clear that such integrals have imaginary part on p 2 , and q 2 - the momentum of 
nucleon and axial current - but there is no imaginary part on p\ - the momentum of 
pentaquark. So we come to the conclusion that such diagrams correspond to the case, 
when there is no G + resonance in the pentaquark current channel (this correspond 
to background of this decay). (Note, that this conclusion don't depend on the fact 
that one ore more of the quark propagators should be replaced by condensate, as we 
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discuss before). The double imaginary part on p\, p\ (i.e. B + resonance and baryon) 
appears only if one take into account hard gluon exchange, and not arbitrary, but only 
those, which connect the quark line, going to axial current vertex with those going to 
an baryon vertex, (so that it provide the moment exchange between these vertexes), 
as on Fig. lb. So we come to conclusion, that if + is a genuine 5-quark state (not, 
say, the NK bound state), then the hard gluon exchange is necessary, what leads to 
additional factor of ot s . We see, that pentaquark width T® ~ a 2 (0|gg|0) 2 , i.e., Te has 
strong parametrical suppression. 

Part 2. 

Let us now discuss the sum rules (7). First of all, of course, the contribution of 
the operator of d = 3 should be accounted. As we discuss in previous section, only 
diagrams like those on Fig. lb give the non- vanishing contribution to double dispersion 
relation. It is necessary to note, that in the sense of discussion before (see [Hj) it is 
quite necessary to keep only those part of these diagrams, which have really imaginary 
part both on p\ and p\ (i.e pentaquark and nucleon 4-momenta square). This mean 
that double Borel transformation (forp 2 andp 2 independently) in (7) is quite necessary. 

In what following, we will denote the result of calculation (i.e. operator d = 3 
contribution to sum rules (7)) as 

R d 3 = l/(X n X e )e^ M - +m2 ^\-l/2)B e B n U QCD (12) 

It is clear, that 

90n = Rd3 + Rd5 + ••• (13) 

and if one confine itself only of d — 3 operator contribution, (as it was done in |15j ) 
then g^ n = R^. Let us shortly discuss the results of |T3j (the higher order operators 
contribution will be discussed in the next section). 

First of all the calculation of the diagrams Fig. lb is technically enough complicated. 
One should pay special attention to extract the terms, which have no imaginary part 
on pentaquark 4-momenta correctly. We perform calculation in the x-representation in 
Euclid space, using standard exponential representation of propagators and the relation 
Be- bp2 = 5{b- 1/M 2 ). 

We neglect in the calculations the effects of s — quark mass, which are very small. 
We also use in calculation the fact, that the ratio Al = M^/Mg, (where M 2 ,M| are 
nucleon and pentaquark Borel masses) should be of order of ratio of the corresponding 
mass square m^/m 2 ,, so we can threats Al as small parameter. But even at this 
simplification the analytical answer is enormous large and is expressed in terms of 
a very large number of different double (and ordinary) integrals and it total size is 
very large (about some hundred terms). That's why, unfortunately, we can not write 
down pure analytical answer, and the last last stage of calculation we prefer to do only 
numerically . (Of course we check, that all integrals converge if Q 2 (= —q 2 ) is not equal 
to zero.) 

Of course, as was discussed before, the result (Rdz) is proportional to quark con- 
densate and strong coupling constant, i.e. a s a, where a = — (27r) 2 (0|gg|0). As a 
characteristic virtuality we chose the Borel mass of nucleon, but one should note, that 
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because a s a 2 do not depend on normalization point, this choice is rather unessential. 
We use the value of 

a s a 2 = 0.23GW 5 , X 2 n = X 2 n * 32vr 4 = 3.2Gev 6 , X 2 = X 2 e * (4vr) 8 = 12Gev 12 (14) 

(see [21], [HI)- To extract the continuum contribution, usually one should write 
down explicitly double dispersion integral, but it is too difficult to do technically in 
our case, so we estimate continuum contribution by a usual factors E — 1 — e s °^ M in 
both pentaquark and baryon channels. The continuum dependence is not strong, we use 
standard value s = 1.5GeV 2 for nucleon and s = 4. — 4.5GeV 2 for pentaquark current 
[TH] . We estimate the inaccuracy of the value of due to continuum contribution 
about 30%. In the paper [T3] only the contribution to the value gf n ( eq.(13)) was 
accounted and the estimation g^ n < lMev was obtained. In this paper we also account 
the contribution of the higher dimension operators we are now going to discuss. 

Part 3. 

There are large number of diagrams, corresponding to higher dimension operators 
contribution (see some examples on figs. 2, 3, 4 for dimension d — 5, 7, 9 correspondingly). 
Numerically the main contribution in our case give diagrams on Fig. 4, (operator of the 
dimension d = 9), proportional to a 3 , where a = — (27r) 2 (0|gg|0)). It is not surprising, 
the same situation appear in many similar cases ( see ^H] for the case of nucleon, or 
[T5] for the case of pentaquark). The reason is clear -each cut quark line lead, from one 
side, to small factor -the quark condensate, but, from the other side, it decreases the 
number of loops, so lead to additional factor of order of 4tt 2 , and this two factors more 
or less compensate each other. That's why one can expect, that main contribution of 
high dimension operator come from diagrams on Fig. 4a-4b. When calculating this 
diagrams, one should also carefully exclude those, which have no double imaginary part 
both on p 2 and p 2 .. Examples of such diagrams which have no double imaginary part 
are shown on Fig. 5 - one can easily check this just in the same way as was discussed 
in the previous chapter (for diagrams on Fig. la). So to estimate dimension d = 9 
operator contribution R d9 one should account only four types of diagrams on Fig.4a,b 
(of course each of them with all possible combinations). Fortunately, in this case it is 
possible to write down analytical answer. If we neglect the small terms, proportional 
to s-quark mass, and also for simplicity suppose strange and light quark condensates 
to be equal, we can write: 

3 so s n 

R d9 = -0.209(A n A e )- 1 e«/ M «+ m «/ JVf ^^- / e'^du J e- s ' M *ds( Pl + p 2 /v) (15) 

o o 

where 

Pl = (-3/8)u5(s)/Q 4 

-1 s 2 - u 2 + Q A - l.hsQ 2 „ u + Q 2 -s 
p2 = W + 2 Qhfl + 7MS ^ 

v 2 = (u + Q 2 — s) 2 + 4Q 2 s, and A n , Ag are defined in (14). 
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In the numerical analysis we suppose again M| = 3M 2 , as for d = 3 operator case 
in previous section. 

On Fig.6a,b the axial constant gf n = R^z + R^, (see eq. (13)), obtained from sum 
rules, is shown for two values of Q 2 : Q 2 = 1.5Gev 2 (Fig. 6a) and Q 2 = 2.5Gev 2 (fig. 6b) 
as a function of Borel mass of nucleon. Thick (upper) line mean the total result for 
gg n , thin line and dashed line - Rd3 and Rdg contributions correspondingly. One can 
see, that Rd9 is smaller than R^ in this region and also the Borel mass behavior of gf n 
is very good, so one can suppose, that sum rule are reliable at this range of Q 2 . 

From the sense of sum rule, it is clear that we can found axial constant gf n only at 
Q 2 not close to zero (about lGev 2 or higher). And really, at Q 2 = lGev 2 Rd9 became 
equal (and even more) than R d3 , so at this value sum rules became very doubtful, and 
at lower values of Q 2 sum rules does not exists. 

On the Fig. 7 the Q 2 dependence of gf n is shown (at Ml = lGev 2 and M| = 3M^ 
-thin (lower) line, and the same forM^ = 1.2Gev 2 - thick (upper) line ). One can see, 
that they are practically identical. 

Really we are interest the value gf n in the limit Q 2 —* 0, which can't be calculated 
directly from S.R., obviously. But one can see from Fig. 7, that Q 2 behavior is found to 
be almost linear so one can extrapolate it to zero. We found averaged (on Borel mass) 
gf n = 0.034 at Q 2 = 0. Of course the accuracy of this value is not high, and highly 
depend on the extrapolation (about a factor 1.5-2 for different more or less reasonable 
extrapolations). Also large inaccuracy appear due to: 

a) the method itself has the accuracy of the order of the SU(3) violation 

c) accuracy in the value of Ag (about 20-30%) 

d) accuracy of sum rule approach, especially for pentaquark case (see, for example, 
discussion in |25]). 

For all this reason, we estimate the value of axial constant gf n (0) with inaccuracy 
about a factor two, so it can varied from 0.017 to 0.068 with central point gf n = 0.034. 
By use of eq (10) one can easily express the pentaquark width in terms of (g^ n ) 2 , and 
we come to conclusion, that our result for T® can vary in the region from 60KeV to 
lMeV. The central value gf n = 0.034 correspond to width about 0.25MeV. Note, 
that our result for width of the pentaquark (with positive parity) don't contradict to 
value 0.75MeV, obtained in also in sum rules, but by quite differ method. 

Our estimation of the pentaquark width also is in the agreement with the result 
jO] (0.36MeV^ with accuracy about 30%), obtained from the ratio between numbers of 
resonant and non-resonant charge exchange events. 

Of course the accuracy of our result is not high, really we predict only the order 
of magnitude of the width, but it seems, that more precise prediction for pentaquark 
width in this approach is impossible. That's why we don't discuss contribution of 
operators d = 5, 7 e.t.c, because they are much smaller than those we have accounted 
and their contribution is within of the our accuracy. 

The main conclusion is, that if the 9 + is genuine states then sum rules really predict 
the very narrow width of pentaquark, less than lMev, (most probably about 0.25 — 
0.5Mev) and the suppression of the width is both parametrical and numerical. 

Author is thankful to B.L.Ioffe for useful discussions and advises. This work was 
supported in part by CRDF grant RUP2-2621-MO-04 and RFFI grant 06-02-16905. 
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Figure 1: examples of the diagrams for d = 3 operator contribution 




Figure 2: Examples of diagrams for d — 5 operator contribution, dashed line mean 
vacuum gluon field, wave line -hard gluon exchange, circles mean derivatives 
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Figure 3: The same as in Fig. 2. for d = 7 operators 
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Figure 4: examples of the diagrams for d — 9 operator contribution 




Figure 6: g$ n dependence on Borel mass for: a) Q 2 = 1.5GeV 2 , b) for Q 2 = 2.5GeV 2 
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Figure 7: Q 2 dependence of gf n ioiM 2 = 1.2GeV 2 -thick (upper) line and for M 2 = 
l.GeV 2 - thin (lower) line 
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